Binding two-loop vacuum-polarization corrections to the bound-electron g factor 
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We commence the evaluation of the one- and two-loop binding corrections to the g factor for an electron 
in a hydrogenlike system of order a^{Za)^ and consider diagrams with closed fermion loops. The one- 
loop vacuum-polarization correction is rederived and confirmed. For the two-loop vacuum-polarization cor- 
rection, due to a specific gauge-invariant set of diagrams with closed fermion loops, we find a correction 
6g = 7.442 (o/tt)^ (Za)^. Based on the numerical trend of the coefficients inferred from the gauge-invariant 
subset, we obtain a numerically large tentative estimate for the complete two-loop binding correction to the g 
factor (sum of self-energy and vacuum polarization). 
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PACS numbers: 31.30.js, 12.20.Ds, 06.20.Jr, 31.15.-p 



I. INTRODUCTION 

The bound-electron g factor has been the subject of intense 
investigations over the past decade, both experimentally as 
well as theoretically. It describes the response of the bound 
electron to an external homogeneous magnetic field and is nat- 
urally different from the g factor of a free electron, due to the 
binding of the electron to the nucleus. Recent measurements 
for hydrogenlike ions with a spinless nucleus in the region 
of low nuclear charge number Z have been reported and dis- 
cussed in Refs. iHllllll. 

For precision experiments with trapped hydrogenlike ions, 
the most important atomic state to be considered is the ground 
state, and we restrict the discussion to the g factor of the elec- 
tronic ground state, for which we write g = 5 (IS). From the 
relativistic (Dirac) theory of the bound electron (which does 
not include radiative coiTections), one obtains 



g = 2-\ [Zaf - \ {Zat + 0{Zaf 
6 6 



(1) 



Here, a is the fine-structure constant, and Z is the nuclear 
charge number. The negative sign of the correction terms of 
higher order in the Za-expansion imply that g < 2 for higher 
nuclear charge numbers Z. Therefore, planned experiments in 
the high-Z region ^ have been termed "2-g"-experiments. 

The quantum electrodynamic (QED) corrections to the 
bound-electron g factor can be expressed as a combined ex- 
pansion in a and Za, where the latter parameter describes the 
strength of the coupling to the nucleus |7]. The first few terms 
in the expansion of the one-loop coiTection to the bound- 
electron g factor (sum of self energy and vacuum polarization) 
in powers of Za read | 



(Za)' 



11.303191 



+ {ZaY 



f ln[(Za)-^] 



(2) 



+ a5o {Zaf + 0{Zaf 



According to commonly accepted conventions, the coefficient 
aso caiTies two indices, the first of which counts the power 
of Za, whereas the second counts the power of the logarithm 

ln[(Za)-2]. 



The two-loop coiTection reads U. 

2 r / (Za)2 



(2) 



-0.656 958 1 
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+ iZaf (3) 



X < flniiZa)-'] 



16.436 842 + hoiZa)^ + 0{Zaf 



Our goal here is to evaluate the contribution to 650 due to a 
subset of the diagrams containing closed fermion loops, and 
to rederive the known result for the vacuum-polarization con- 
tribution to flso. We recall that according to Fig. 21 of Ref. llj], 
the number of two-loop diagrams contributing to the g factor 
is large, and the particular diagrams considered here form one 
of the most straightforward gauge-invariant subsets in two- 
loop order. As the whole set of two-loop diagrams can be 
broken up into smaller gauge-invariant subsets, the evalua- 
tion could be initiated by considering gauge-invariant subsets. 
Since the diagrams for the g factor are related to those for the 
Lamb shift (except for the additional presence of an external 
magnetic field), a byproduct of our calculations is a confir- 
mation of results obtained previously for the contribution of 
corresponding diagrams to the two-loop, binding correction 
to the Lamb shift 1 1 1] of order a^ {Zaf. 

This brief communication is organized as follows. After 
a discussion of the one-loop correction in the order a [Zaf 
(see Secinii, we describe the two-loop calculations in the or- 
der {Zcxf' in Sec.|III] Conclusions are drawn in Sec. lIVI 



II. ONE-LOOP CORRECTION 

First, we would like to rederive the leading vacuum- 
polarization coiTection to the bound-electron g factor of order 
a(Zaf. To this end, we recall that for the interaction of an 
electron with a constant magnetic field, one can derive the 
following, effective Hamiltonian based on long-wavelength 
quantum electrodynamics [12] for the interaction of an elec- 
tron with an external static magnetic field B, 



1 



2m 4™-^ 



-^{r-^V) \ , (4) 
12 m-^ / 



where p is the bound-electron momentum, m is the electron 
mass, and V is the total static potential felt by the electron. 
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This potential can be either the Coulomb potential, which we 
denote by Vc in the following, or a vacuum-polarization cor- 
rection SV. 

We now briefly recall how to evaluate of the one- 
loop vacuum-polarization correction based on the effective 
Hamiltonian ^ and on well-known formulas for vacuum- 
polarization effects. Indeed, we use the well-known Uehling 
approximation for the vacuum-polarization potential and 
identify the potential in (|4|i as — > SV Vjj, 



TT 



dv 



exp (—A r) 



-Za 



with A = 2 m/Vl - v^. 

The first correction Ei to the spin-dependent magnetic-field 
interaction energy (and thus to the g factor) is obtained if we 
replace V ^ 1% in the third term in brackets in Eq. (|4|i. 
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(6) 



where \<f>) denotes the nonrelativistic atomic ket vector cor- 
responding to the atomic state under investigation (here, the 
ground state). Of course, the rightmost term in Eq. (|6]l is eval- 
uated on the bound-state wave function, but we write it as 

being proportional to (^a ■ B^, where it is understood that for 

a S state, the spin is either pointing up or down. This means 
that the expectation value on the right-hand side is to be evalu- 
ated using the spin degrees of freedom only, and it is therefore 
denoted by a simple bracket. Because Ei is a first-order spin- 
dependent energy correction in a uniform external magnetic 
field, it can be related directly to a correction to the g factor 
For this purpose, we write the interactions as multiplicative 
corrections to the normalized interaction — -r^ <j ■ B; the lat- 
ter leads to a g factor of unity. 

The correction Ei, which is a first-order correction, now 
has to be supplemented by some second-order effects. Let us 
therefore consider the case where V in Eq. ^ represents the 
Coulomb potential Vc- In order to evaluate the second-order 
effects, we investigate the Uehling correction in conjunction 
with the second and the third term in round brackets in Eq. (01, 
which represent corrections to the a ■ B interaction of relative 
order {Za)'^. The perturbation to the wave function induced 
by the leading-order interaction — etr • B/{2m) = —gea ■ 
B /{Am) vanishes. 

The first of the nonvanishing second-order effects is ob- 
tained by considering a second-order perturbation involving 
the Uehling potential and the second term in brackets in 



E-2 
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The second of these is obtained by considering again the third 
term in brackets in (|4|i, but this time acting on the Coulomb 
potential V in second-order perturbation theory. 



E?, 



--2(6 — 




[{r-V)Vc]S-B 
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(5) Taking into account the Hellmann-Feynman theorem. 



(8) 



(9) 



the sum of the corrections Ei + E2 + E^ leads to the known 
result lH 
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(10) 



In the lowest-order in the Za expansion, Eq. (|5]l then imme- 
diately leads to the leading-order vacuum-polarization correc- 
tion to the g factor [9], 



Sg^^ = l{Zaf[-^^ 



(11) 



where the index U reminds us of the Uehling potential. 

We now consider the wave function slope and the a {Za)^ 
correction. According to Schwinger's textbook ifisll . one can 
obtain the vacuum-polarization correction of order a (Za)^ 
to the Lamb shift by considering the slope of the bound-state 
wave function at the origin. This holds equally well for the 
g factor The reason is that the bound-state wave function 
decays exponentially as exp{~Zamr) whereas the Uehling 
potential decays much faster, namely according to Eq. (|5]) as 
cxp(— Ar) where A is of the order of the electron rest mass. In 
the resulting product 

Vu{r) ~ exp(— ZaTOr — Ar) 

= exp(-Ar) (1 - Zamr + 0{r^)) , (12) 

one can thus expand in the first argument of the exponential, 
using A ^ Zarri- The correction term 1 — Zamr corresponds 
to the slope of the wave function at the origin. A straightfor- 
ward evaluation gives the following vacuum-polarization cor- 
rection for the ground state: 



Sgvp 



-15 



{Za) 



(13) 



which includes the correction of relative order Za. We here 
confirm the result of Ref. [14]. For completeness, it is use- 
ful to recall the corresponding one-loop vacuum-polarization 
correction to the Lamb shift, which reads jisll 

5E^^l = ^ {Za)'m (^-1 + %{Za)^ , (14) 
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This concludes the consideration of the one-loop vacuum- 
polarization correction of order a (Za)^ to the g factor 




FIG. 1: Feynman diagrams for the two-loop vacuum-polarization 
corrections to the bound-electron g factor. The first of these is the 
loop-after-loop Uehling vacuum-polarization correction and gives 
a contribution of — 3687r/1701 in units of (a/yr)^ (Za)^ for the 
g factor. The sum of the second and third diagram (upper row), 
which are iterated one-loop perturbations, yields a contribution of 
— 85l7r/6804 in the same units. The two last two diagrams (lower 
row) are Kallen-Sabry diagrams. They lead to a <; factor correction 
of TT (125176/19845 + 832 ln2/189 ~ 4007r/189). 



III. TWO-LOOP CORRECTION 

We now discuss the two-loop calculation. The leading- 
order result 1 10] follows as we replace the Uehling potential 
Vu in Eq. ^ by the KaUen-Sabry term. It reads 



JKS 



(15) 



The Kallen-Sabry term is a genuine two-loop effect, and one 
might wonder why the iterated one-loop diagrams do not also 
contribute in the order a^{Za)'^ (these correspond to the sec- 
ond and third diagrams in Fig. [1]). However, the first-order 
correction to the wave function induced by the leading-order 
magnetic interaction —ea ■ B/{2m) vanishes, and the re- 
maining contribution due to the iterated Uehling term can- 
cels explicitly in third-order perturbation theory, because the 
term with the magnetic interaction "in the middle" cancels 
against the derivative term obtained by considering the deriva- 
tive of the second-order Uehling correction with respect to the 
bound-state energy. 

For some of the diagrams in Fig. [T] the g factor correc- 
tion of order {Za)^ can be obtained by expanding the wave 
function about the origin, as it was done for the one-loop the- 
ory. An example is the first diagram in Fig. [T] which can 
be expressed as the expectation value of the loop-after-loop 
Uehling potential, evaluated with wave functions perturbed 
by the magnetic interaction. For the other diagrams, the cal- 
culation is more complicated. In particular, since the iterated 
Uehling correction (second and third diagram in the first row 
in Fig. [Til contributes at the order of a^{Za)^, one cannot 
avoid the complete calculation of the (first-order) perturba- 
tion to the wave function by the Uehling correction, which 



involves exponentials, exponential integrals, logarithms and 
powers of the radial variable. E.g., the magnetic interaction in 
the middle vertex demands a further integration over the radial 
coordinate. The further calculation proceeds along the lines 
outlined in Ref. 1, 1 1,1 for the two-loop vacuum-polarization 
corrections to the Lamb shift. 

We finally obtain for the two-loop binding contribution 
Sgyp due to the diagrams in Fig.[T] 



— + {Za)7r 



/ 1420807 
V 238140 



832, „ 400 

In 2 TT 

189 189 



= [Zaf [-4.049 + 7.442 (Za)] . (16) 

The numerical coefficient of the (ZQ;)-correction is rather 
large, mainly because it has a factor tt in the numerator 

Just as for the one-loop calculation, it is useful to compare 
our results to those for the Lamb shift, selecting the corre- 
sponding set of diagrams. For the Lamb shift, we can identify 
the diagrams corresponding to those in Fig.[T]by simply elim- 
inating the interaction with the external magnetic field. The 
resulting diagrams after this removal operation are equivalent 
to the diagrams labeled as IV and VI in Re f. 11 hi- The corre- 
sponding contribution to the Lamb-shift is iflUllql 



52 25 

m 2 TT 

63 63 



82 , , /7421 

h {Za)TT 

81 ^ ' V6615 



a 



(Zaf [-1.012 + 1.405 {Za)\ , (17) 



and we have verified it using our approach. This concludes 
our two-loop vacuum-polarization calculations. 



IV. SUMMARY 

In this brief report, we describe the evaluation of a part 
of the binding, vacuum-polarization correction to the bound- 
electron g factor. The vacuum-polarization corrections repre- 
sent a preparatory calculation for the self-energy corrections, 
which are much more difficult to evaluate. In view of the mul- 
titude of terms generated in comparison to the corresponding 
self-energy correction to the Lamb shift of order a{Za)^, and 
in view of the additional complexity of the calculation due to 
the added external magnetic field, considerable difficulties are 
expected. 

It may, akeady at this point, be permitted to speculate a 
little about the magnitude of the complete correction to the 
g factor of order a'^{ZaY\ which is less than an estimate 
but perhaps more than just guesswork. Namely, we observe 



there appears to be a rather universal factor in the range of 
3.5 ... 5.5 by which the g factor coefficients of a given or- 
der in the Za-expansion are larger than the corresponding 
Lamb shift coefficients for the ground state. In particular, 
we compare in the order a(Za)'^, the coefficient —16/15 in 
(fTsT l to the coefficient —4/15 in Eq. (fT4l l (the g factor coef- 
ficient is larger than the Lamb shift coefficient by a relative 
factor four). At relative order Za, the relative factor is 5.33 
(the coefficients are 5tt/9 versus 57r/48). For the g factor at 
two-loop order, the relative factor at a'^{Za)^ is 5.30, as evi- 
dent from Eqs. ( fT6] l and ( fTTj i. A factor in the range 3.5 ... 5.5 
also appears for the self-energy corrections. We recall that the 
complete two-loop correction to the Lamb shift in the order 

SE^ ^ -21.55 {Zoifm. (18) 



4 

Our "educated guess" for the complete correction to the g fac- 
tor thus is 

^g^^) = C {Zaf , -118<C<-75. (19) 

The magnitude of this estimate of the coefficient generates 
obvious interest. 
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